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Abstract 

The gravitational analogue of the electromagnetic Meissner effect is 
investigated. Starting from the post-Newtonian approximation to general 
relativity we arrive at gravitational London equations, predicting a grav- 
itational Meissner effect. Applied to neutron stars we arrive at a London 
penetration depth of 12km, which is about the size of a neutron star. 

PACS number(s): 04.60, 74.25.H, 97.60. J, 04.70 



I Introduction 

There is strong evidence that neutron stars are superfiuid and maybe super- 
conducting indicating that macroscopic quantum effects play a significant 
role in neutron stars. In addition neutron stars are very strongly gravitating 
objects. We therefore would like to investigate the possibility of macroscopic 
quantum gravity effect in the form of a gravitational Meissner effect. 

In the post-Newtonian approximation Einstein's theory of general relativ- 
ity, gravity looks like a U(l) theory and is very similar to the electromagnetic 
Maxwell theory || . This allows us to follow closely the approach of London [|| |j 
in deriving gravitational London equations, predicting a gravitational Meissner 
effect. Not knowing the underlying microscopic theory of quantum gravity we 
have to resort to this semi-classical approach. We find the gravitational London 
penetration depth to be about 12 km for a typical neutron star of radius 15 
km. Furthermore, as the density of the neutron star increases, the London pen- 
etration depth decreases, suggesting that quantum gravity aspects may become 
significant as the neutron star approaches collapse. 

Macroscopic quantum aspects are best described by a BCS type theory Q. 
However, since BCS needs a quantum theory of the underlying interaction, and 
since we do not have a microscopic theory of quantum gravity at this moment, 
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a gravitational BCS type theory is not within our reach. Despite this, we can 
certainly make a plausible argument for what happens at the microscopic scale. 
Inspired from the Cooper-pair forming electrons in the BCS theory of supercon- 
ductivity, the Interacting Boson Model (IBM) M, fj| assumes nucleons to form 
pairs inside nuclei. These pairs are bosons and by only looking at the interaction 
of those bosons among themselves, one can explain a wealth of nuclear data to 
a high degree of accuracy. Neutron stars consist mostly of neutrons and some 
protons and for low enough temperatures (< 10 9 K) they are expected to pair 
up. Experimental evidence for the pairing comes from the theory of cooling of 
neutron stars, where the reduction of the specific heat due to extensive pairing is 
essential § . In addition, bose-condensation can occur and superfluidity as well 
as superconductivity are consequences. Evidence for superfluidity can be found 
in the glitches in the timing history of pulsars |jl], [|. Therefore, although we do 
not have a microscopic theory of quantum gravity, we have enough evidence jus- 
tifying a phenomenological investigation of macroscopic quantum gravity effects 
following London's approach. 

Symmetry breaking is a very convenient way to look at macroscopic quantum 
effects. To break a symmetry one needs a bose-condensate. What causes the 
bose-condensation and which symmetry is broken are two entirely unrelated 
problems. Any mechanism leading to an attraction between the fermions at 
the Fermi surface is equally well suited for producing the bose-condensate. For 
example, in the standard BCS theory of superconductivity, the attraction is 
caused by the phonon interaction, the symmetry broken is the electromagnetic 



U(l) symmetry 10 1. In the theory of superfluid helium 3 the attraction is 
caused by van der Waals forces, the symmetries broken are translational p"2] ] and 
rotational in spin and orbital space as well as gauge symmetry. In the theory 
of the Interacting Boson Model, the attraction is caused by the strong nuclear 
force, the symmetry broken is related to angular momentum and Runge-Lenz 
vector Q. In the theory of superconductivity of protons inside neutron stars, 
the attraction is caused by the strong nuclear force, the symmetry broken is the 
electromagnetic U(l) symmetry. Therefore, we see no reason why to exclude 
the scenario, where the attraction is caused by the strong nuclear force and the 
symmetry broken is gravitational. 

This work is also an attempt to probe certain aspects of quantum gravity. 
The microscopic approach to quantum gravity has been poised with serious 
problems, and no satisfactory theory has emerged as of yet. The macroscopic 
approach seems much more tractable, and maybe following the theory of su- 
perconductivity from London to Ginzburg and Landau to Bardeen, Cooper and 
Schrieffer will give us another route to learn more about quantum gravity. 

II Gravitational Maxwell Equations 

Let us review the post-Newtonian approximation fl3], nil. For weak gravi- 
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tational fields and low velocities, i.e., 



U/c 2 ~ («/c) 4 ~ (p/pc 2 ) 2 ~ (n/c 2 ) 2 « 1, 

Einstein's field equation can be written in a form resembling Maxwell's equations 
of electrodynamics. We are mostly interested in magnetic-type gravity, and 
therefore we shall use the truncated and rewritten version of the parametrized- 
post-Newtonian (PPN) formalism given by Braginsky et. al. JL5| . Scalar and 
vector potentials are defined as 

= -((7 + 2*) and Aj = -Ivj - ^Wj, (1) 

where {/, \&, V}, and Wj are defined in chapter 39 of Misner, Thorne and Wheeler 
[ p"4| . This allows us to define a 'gravitoelectric' field g and a 'gravitomagnetic' 
field H through 

1 <9A 

g = -V<f> and H = V x A. (2) 

c at 

In this weak field, slow motion expansion of general relativity, g contains 
mostly first order corrections to flat space-time, and H contains second order 
corrections. For the solar system, g is just the normal Newtonian gravitational 
acceleration, whereas H is related to angular momentum interactions and effects 
due to H are about 10 12 times smaller than those due to g. 

In the above approximation, Einstein's field equations can be rewritten in a 
form very much resembling that of Maxwell's equations 

1 r)l-T 

V • g = -4nGp a + O (c- 2 ) , Vxg = --— (3) 

' c at 

V-H = 0, VxH = 4f-47rG^ + J|j, (4) 

where po is the density of rest mass in the local rest frame of the matter and v 
is the ordinary velocity of the rest mass relative to the PPN coordinate frame. 
The specific internal energy II and the pressure p have been neglected since they 
only contribute O (c~ 2 ) effects to the electric field. 

The equation of motion of an uncharged particle is identical to the electro- 
magnetic Lorentz force law [|[ |l5| , 

^g + JvxH + 0(c- 2 ) (5) 
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Ill Gravitational London Equations 



With these gravitational Maxwell equations we can derive the equivalent of 
London's equations for gravity. The derivation is completely analogous to the 
electromagnetic case j|, ||, 16, |l7|. The Lorentz force law tells us that given a 



gravitational electric field there will be a resultant acceleration. Multiplying by 
po, the matter density, and defining the matter current j = pov, we arrive at 

dt = Pag - (6) 
This is sometimes referred to as the first London equation. Now substitute this 
into the right equation of (|jj|), the equivalent of Faraday's law, we obtain 

One possible solution to the above equation is obtained by integrating and 
setting the integration constant to zero. This leads to the second (and most 
famous) London equation 

Vxj = -^H, (8) 

which predicts the Meissner-Ochsenfeld effect ]l8j . 

In the Ginzburg-Landau theory of superconductivity, it is shown that setting 
the integration constant in equation ((7) to zero is equivalent to the assump- 
tion that the significant spatial variation of the order parameter ip = \ip\ e l< ^ is 
through the phase (j), and not the magnitude \tp\ For the density of 

Cooper pairs this means that it cannot vary appreciably from its uniform ther- 
mal equilibrium value. For a superfluid neutron star in thermal equilibrium, 
where the pairs can flow, but not accumulate or be destroyed, this should be 
the case. 

To obtain the London penetration depth, we take the static part of the right 
equation of (||) 

„ TT 16vrG . 

VxH = j (9) 

c 

and take the curl of it. Using the second London equation (^) we arrive at 

V 2 H=16^^H, (10) 
where we immediately identify the gravitational London penetration depth 

Ai = yy 1,2 = 5177 x 1012 p " 1 '"- (u) 

where the density po is measured in kg/m 3 and is given in meters. 
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IV Coherence Length 



For our derivation of the London equations to be valid we must require that 
the London penetration depth is much larger than the coherence length £o 

E- 

A L »£o, (12) 
where the coherence length is defined through 

« 13 » 

with vf the Fermi velocity and 2A the pairing energy or energy gap. We are 
only interested in a rough estimate of the upper limit of £ - An upper limit for 
the Fermi velocity obviously is the speed of light. If the pairing is caused by the 
strong nuclear force, as we assumed, then a lower limit on A is a few hundred 
thousand keV With these numbers we find for the coherence length 

Co < lCT 13 m, (14) 

which is larger than the separation between two nuclcons, but much smaller 
than the London penetration depth. Therefore, our derivation of the London 
equations was justified in the given approximation. 



V Discussion 

We notice that for regular densities occurring in the universe the gravita- 
tional London penetration depth will be a very large number. However, inter- 
esting enough, for neutron stars with densities of about p^s — 2 x 10 17 kg/m 3 , 
we obtain a London length of 12km, which is slightly smaller than the radius of 
a neutron star with that density. The London length is inversely proportional 
to the density, meaning with increasing density the London penetration depth 
decreases. 

For neutron stars this means that the onset of a gravitational Meissncr effect 
is predicted. In essence it implies that the gravitational magnetic field caused 
by the huge angular momentum of the neutron star, will be expelled from the 
center of the neutron star. This could be accomplished through induced matter- 
supercurrents in the outer layers of the neutron stars, creating counter-magnetic 
fields to expel the gravitational magnetic field from its interior, in analogy to 
the electromagnetic case. 

A more detailed description is not possible in the approximation used, since 
the post-Newtonian approximation is only a rough approximation if used for 
neutron stars. We are currently working on the second order post- Newtonian 
approximation and numerical simulations |2l| . 
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VI Conclusion 



More significant than the expulsion of the gravitational magnetic field, how- 
ever, may be the fact that the Meissner effect is a macroscopic quantum effect, 
and in essence, the gravitational Meissner effect is a macroscopic quantum grav- 
ity effect. There is experimental evidence that neutron stars are superfluid as 
well as superconducting. This indicates that macroscopic quantum effects do 
play a significant role for neutron stars. Furthermore, neutron stars are the 
strongest gravitational objects observed with certainty. Under the assumptions 
used here, we find that macroscopic quantum gravity effects may play a signifi- 
cant role for neutron stars. 

Since the work of Oppenheimer and Volkoff P2| , we know that neutron 
stars will be stable only up to a certain point. After that the gravitational 
pressure overcomes the pressure of the neutron Fermi gas and continued collapse 
is predicted. This leads to increasing densities, and in our model to decreasing 
London penetration lengths, indicating that gravitational quantum effects might 
become more and more important. Therefore, we would like to emphasize that 
the prediction of Oppenheimer and Volkoff is classical in nature, and our work 
indicates, that a quantum gravitational treatment may lead to new results. 
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